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We study the nonlinear feedback in a fermion-boson system using an extension of dynamical mean- 
field theory and the quantum Monte Carlo method. In the perturbative regimes (weak-coupling and 
atomic limits) the effective interaction among fermions increases as the width of the boson dispersion 
increases. In the strong coupling regime away from the anti-adiabatic limit, the effective interaction 
decreases as we increase the width of the boson dispersion. This behavior is closely related with 
complete softening of the boson field. We elucidate the parameters that control this nonperturbative 
region where fluctuations of the dispersive bosons enhance the derealization of fermions. 

PACS numbers: 71.38.+i, 63.20.Dj, 63.20.Kr, 71.30. +h 



I. INTRODUCTION 

Interacting fermion-boson systems are very impor- 
tant in condensed matter physics and have been studied 
intensively.!!! They are directly relevant to the description 
of electron-lattice interaction. Other problems can be 
mapped onto interacting fermions and bosons hx means 
of the Hubbard-Stratonovich transformationotl While 
the problem of a single fermion interacting with a bet 
son field, i.e., the polaron problem, is well understood, El 
a lot less is known about the many-fcrmion problem in in- 
teraction with a boson field; it is a full interacting many- 
body problem which is only tractable analytically in the 
adiabatica and the atomic limits. ultllil 

In this paper we revisit the interacting and dis- 
persive fermion-boson problem using dynamical mean- 
field (DMT) theory.Q This method reduces the quantum 
many-body problem to a quantum impurity model obey- 
ing a self-consistency condition. This method has been 
useful in describing strong coupling problems such as the 
Mott transition. There are several motivations for our 
work. 

First, a DMF treatment of the bosonic and fermionic 
degrees of freedom taking into account the boson disper- 
sion, requires an extension of the DMF equations where 
the bosonic propagator degrees of freedom are deter- 
mined self-consistently. This represents a new type of 
self-consistent DMF equation, which so far has not been 
investigated. These equations are relevant to many prob- 
lems, electron-phonou. interactions, fermions interacting 
with spin fluctuations!] or among themselves via the long- 
rangecl-Coulomb interactions^] and to the boson-fermion 
model.EJ 

Second, while the Mott transition in the Hubbard 
model is well understood using DMF methods, it is inter- 
esting to understand how it is modified by the variation of 
the frequency of the mode that mediates the interaction, 
or how the results are changed by the electron-phonon 
interactions. The approach discussed in this paper is a 
first step in this direction. 



Finally, phonon dispersion effects are relevant to many 
systems. The Jahn- Teller or breathing-type phonons, 
for instance, seen in manganese oxides should be disper- 
sive due to intersite coupling. A distortion of a MnCv 
octahedron affects distortions of the neighbor octahe- 
dra, since the Mn06-octahedra share their oxygen atoms 
which leads to an intersite coupling. This may be rel- 
evant to fascinating-, oederings of lattice and charge in 
doped manganites.tirfc-3 

We study the mutual feedback of fermionic and bosonic 
degrees of freedom in a very simple system of fermions in- 
teracting with one branch of bosons at half filling. How- 
ever the methodology can be extended to other prob- 
lems where similar DMF equations occur such as electron 
problems with long-ranged Coulomb interactions and the 
competition of magnetic order and the heavy fermion 
state, and to the boson-fermion mixture of high temper- 
ature superconductivity. 

This paper is organized as follows. In the next section 
we discuss how DMF theory needs to be extended to 
fully include the feedback effects through fermion-boson 
interaction. Quantum Monte Carlo (QMC) method is 
introduced to solve the DMF equations for all param- 
eter regions. We also discuss some technical points of 
the QMC relevant to this problem. The formalism is ap- 
plied to demonstrate effects of boson dispersion in a wide 
region of parameters and the results are summarized in 
Sec. III. In Sec. IV, we discuss our main result: The ex- 
istence of two distinct regimes of the DMF solutions. In 
the first regime, the feedback effects increase the fermion- 
boson coupling. In the second regime, which is strongly 
fluctuating, the boson dispersion accelerates the dereal- 
ization of fermions. Complete softening characterizes the 
crossover between these regimes. Section |v| is devoted to 
summary. 
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II. DYNAMICAL MEAN-FIELD FORMALISM 
AND HAMILTONIAN 

In this work, we discuss such feedback effects caused by 
the fermion-boson interaction using DMF theory. DMF 
theory provides a local view of a many-body problem 
in terms of an impurity model which satisfies a self- 
consistency condition.0 For general fermion-boson prob- 
lems with a local interaction, the local action has the 
form 

ScS = [ drdr' £ cl(r)g^( T - r')c a {r') 

+ f dTdT'Y, X AT)V^{T-T')x v {T') 
J v 

+ / dr ^ K ia2 uC f ai {T)c a2 (T)x u (T), (1) 
aia2f 

where Go and £>o are the bare impurity Green's func- 
tions for fermion and boson, respectively, which contain 
the dynamical information of the integrated other sites. 
Here c a is the fermion annihilation operator and x v is 
the boson field. X aia2l/ denotes the coupling between 
fermions and bosons. The index a (v) denotes internal 
degrees of freedom of fermions (bosons) such as spins or 
orbitals of electrons (normal modes of phonons). We do 
not explicitly write the contribution from fermion inter- 
actions such as the Coulomb interaction since we focus 
on the effects of boson dispersions in this paper. How- 
ever the action ([!]) is quite general which contains such 
fermion interactions through the Hubbard-Stratonovich 
transformatiorolj with continuous fields. Of course, al- 
ternatively, one can include additionally the fermion in- 
teractions according to the DMF theory for the Hubbard- 
type modelsE 

The full Green's functions are related to the bare ones 

by 



G a X {iU n ) = G a( iuj n) ~ Sq(«W„), 



(2) 
(3) 



at each Matsubara frequency oj n = (2n + for 
fermion and lo„ = 2mr/P for boson, respectively (n is 
an integer). (3 is the inverse temperature. £ and n are 
the self-energy for fermion and boson, respectively. The 
Green's functions for both fermion and boson are deter- 
mined in a self-consistent way. This is achieved by the 
following set of self-consistency conditions, 

Qa = ^ \ib} n + /i — eg a — E a (iu> n )] , (4) 

q 

T> u = y^ j [(«w„) 2 - Jq u - U u (iLu n ) , (5) 



where eq a and ujq v give the dispersion relations for 
fermions and bosons as a function of the wave number 



q, respectively, fi is the chemical potential to control 
the density of fermions. Here the bosons are described 
as harmonic oscillators. The condition (|^) is modified 
according to the property of boson degrees of freedom. 

Previous studies of other models have indicated that 
the results of the DMF theory can give useful insights into 
three-dimensional systemsE We have therefore taken the 
dispersions eg and LJq that correspond to a semicircular 



density of states, see the details in Sec. Ill A. These 
DMF equations are exact for a model where the fermions 
and the bosons have random hopping on lattice sites. 

The self-consistency loop is closed as follows: The ef- 
fective action ([!]) is solved for given bare impurity Green's 
functions Go and Vq to obtain the full Green's functions 
G and V. The self-energy E and n are calculated by 
the relations (Q) and (^) , and used to obtain the Green's 
functions through the self-consistency conditions (^) and 
(||). New bare impurity Green's functions are calculated 
by the relations (||) and (||) again. This loop is iter- 
ated until all the quantities are converged. In this way, 
both fermionic and bosonic dispersions are renormalized 
through the fermion-boson interaction, and the mutual 
feedback effects are fully included. 

The above DMF equations assume that no symmetry 
breaking is present in the system although the extension 
to phases with broken symmetry is straightforward. And 
they can be derived from an electron-phonon model: 



Ti. = Hf + H.b + Hi, 



where 



a q 



(6) 



(7) 



(8) 



v q 



fti — E E M a i^ c p+q ai c pa-2 [a-qv + a -qi/j > ( 9 ) 

where aq is the boson annihilation operator which is re- 
lated with the boson field by xq — (2Muiq)^ 1 ^ 2 (aq + 
al q ). 

The model (|^) has been intensively studied using DMF 
methods in the iiwft of zero boson dispersion, i.e., in the 
Holstein modelHEl Bosons with the same index v have 
a same frequency (Einstein phonons) and the fermion- 
boson coupling is local as 

W b = w „ Y \ a \v a q V + o ) 
v q \ ' 

= E^£ (4, (io) 

v i 
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where the index i denotes a lattice site. In the ground 
state, the possibility of charge-ordered or supercon- 
ducting-, states has been intensively discussed for this 
model.Ejliia Above the critical temperatures of these 
states, the crossover behavior is observed from the Fermi 
liquid with a mass enhancement in the weak coupling 
region to the so-called polaron which is a combined ob- 
ject between, fermion and boson in the strong coupling 
region.™!^ 

It is instructive to compare the present framework with 
the DMF theory for the problem without the boson dis- 
persions such as the Holstein model. If bosons have no 
dispersion, that is, all ujq take the same value u>o inde- 
pendent of q, Eq. (JsJ) is rewritten as 



V = [(iu)„ 



(12) 



Although the full Green's function V contains a feed- 
back effect in the self-energy II, the bare impurity Green's 
function 2? is fixed at the noninteracting Green's func- 
tion given by 



V, 



free 



(13) 



throughout the self-consistency iterations when we start 
from the solution Vq — Vq cq . This is equivalent to the 
ordinary DMF theory for tbe-Holstein model which does 
not need Eqs. (||) and (|])£§E2l Compared to this, for the 
cases with finite bosonic dispersions, the bare impurity 
Green's functions Vq is renormalized from V f Q TCC in the 
iterations in our formalism. 

The renormalization of Vq plays a crucial role be- 
cause Vq is related to the effective interaction between 
fermions. If we integrate out the boson variables x in the 
Hamiltonian, the effective interaction between fermions 
takes the form 

]T J dTdT'cl(T)c a2 (T)UMr-ryi 3 (T')c ai (r'), 

(14) 

where 



with the width At (At = (3/L). Continuous variables 
Xvi = x v {t{) (t; = IAt, I = 1,2, • • -,L) are randomly 
updated to x' vl with the probability 



n 



det Q a exp \~AtB{x' v1 )\ 
detg' a exp[~ATB(x vl )Y 



(16) 



where B(x u i) = J2j=i x ^j' D ouji x ^i witn v ov 3 i = ^Ou(rj 
ti). The fermion Green's functions Q and Q' are calcu- 
lated by the standard algorithrrcil for the configurations 
with x v \ and x' vl , respectively. 

In actual QMC samplings, we consider both local and 
global updates for the continuous fields x v \. The local 
update consists of sequential updates of the fields on each 
discretized point; a change from x v \ to x' v , = x v i + rS 
is attempted where r is a random number between — 1 
and 1 and S is a given amplitude. The global one is a 
simultaneous movement of all the fields at a same amount 
r5. The latter becomes important especially in the strong 
coupling region and/or at low temperatures where the 
fields x show some orderings or are nearly ordered. The 
update amplitude S is chosen to give an appropriate value 
of the acceptance ratio which is defined as the ratio of 
the number of accepted samples to the total number of 
trials. 

The QMC calculations generally have the negative sign 
problem; the MC weight jl^ ) can be negative for the 
general action (|l|), which leads to numerical instability 
in the QMC measurements. However, if fermions couple 
to bosons only in the diagonal form, that is, the coupling 
parameter A QlQ2t , is nonzero only for the case of-ai = a?2, 
the MC weight (|l6|) becomes positive definite.E3 In this 
case, there is no negative sign problem for all parameters. 

There are two sources of errors in the QMC calcula- 
tions. One is a systematic error due to the discretization 
of the imaginary time, and the other is a statistical error 
from the random sampling. The former error is known to 
be proportional to (At 2 ). Measurement is divided into 
several blocks to estimate the latter statistical error by 
the variance among the blocks. The size of each error 
depends on a specific form of models and parameters. 



U eS {T) = X 2 V {t). 



(15) 



In the absence of the boson dispersion, since Vq is un- 
changed through the self-consistency loop as mentioned 
above, the effective interaction ( |l5| ) is also unrenormal- 
ized. On the other hand, Vq is renormalized in our for- 
malism for finite dispersions, which means that the effec- 
tive interaction between fermions is renormalized by the 
mutual feedback of the fermion-boson coupling. 

There are several techniques to solve the effective im- 
purity problem with the. action ([!]). In this work, we 
employ QMC methodErEil because it is an unbiased cal- 
culation and suitable to investigate all the parameter 
regions beyond perturbative regimes. In the QMC ap- 
proach, the imaginary time is discretized into L slices 



III. RESULTS 

A. Model and Parameters 

We apply the new DMF framework proposed in the 
previous section to a case that the general Hamiltonian 
(^) contains two species of fermions and one branch of 
bosons. We set the mass M = 1. The model is a straight- 
forward extension of the Holstein model to include dis- 
persive bosons, whose fermion-boson interaction is given 
by 



Hi - -A 



(17) 
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where the index a takes two values like spin degrees of 
freedom of electrons. The interaction is diagonal in the 
fermion index a so that the QMC does not suffer from the 
negative sign problem as mentioned in Sec. |l[ We take 
the coupling parameter A to be positive, which favors a 
doubly-occupied or an empty state on each site. Note 
that the model has the particle- hole symmetry at fi = 0. 

The boson dispersion is taken into account through Eq. 
(||) in the present framework. We replace the summations 
over the wave number q in Eqs. (|J) and ([5]) by the energy 
integrations as 



Q{iu) n ) = 



V(iuj n ) 



Dp(e) de 



■H — e — T,(iu n ) ' 
D^{e) de 



liilUJn) ' 



(18) 
(19) 



where D-p and Z?b are the the density of states for fermion 
and boson, respectively. In the following calculations, we 
assume semicircular density of states as 



D F (e) 



ttW- 



■VW 2 



e 2 



(20) 
(21) 



where W is the half-bandwidth of the fermion density of 
states which is taken as unity hereafter (W = 1); ojq and 
toi are the center and the half-bandwidth of the boson 
density of states, respectively (iv > 0, uj — u)i > 0). For 
the semicircular density of states, the integrations fl8| ) 
and (H) are performed analytically^ which give 



Q = 



( - VC 2 - 4t 2 



2t 2 



T> = ~ 



(22) 
(23) 



where £ = iui n + /i and £± = £,±u>o with £ 2 = {iuj n ) 2 — n. 

The shape of the boson density of states near the bot- 
tom is important because bosons at the band edge can be 
easily excited and strongly interact with fermions. The 
semicircular density of states ( |2~l| ) has an e 1 / 2 -singularity 
which is expected for bosons with ordinary cosine disper- 
sions in three dimensions. Therefore we believe that the 
following results are qualitatively unchanged in realistic 
three-dimensional models. Results would be different for 
the two-dimensional density of states which has a step 
discontinuity at the band edges and results in very dif- 
ferent DMF solutions. 

In the absence of the boson dispersion {u)\ = 0), 
the model with the interaction ([l7]) (the ordinary Hol- 
stein model) shows a charge ordering around half-filling 
(p = 0) and superconductivity in doped regions at very 
low temperatures.EjtJ In the following, we examine ef- 
fects of boson dispersions in the low temperature region 



above and around these transition temperatures at half- 
filling (/i = 0) assuming no symmetry breaking. The 
calculations are mainly performed at j3 = 8. We take 
At = 1/4 for which all the measured quantities are con- 
verged to the limit of At — > within the statistical 
errors. We have typically run 1,000,000 MC steps for 
measurements; one MC sampling means a set of a sweep 
of local updates over the whole discretized points and a 
global update. Convergence in the self-consistency loop 
is usually rapid; typically 10 iterations are required to 
converge within the statistical errorbars when we start 
from the noninteracting Green's functions. However in 
the strong coupling case, the iteration often suffers from 
an oscillation between two solutions. To avoid the os- 
cillation, we make the iteration proceed by mixing the 
previous solutions. 



B. Dispersionless Boson 

First, we reconsider the limit without the boson disper- 
sion, that is, uj\ = 0. In this case, we use the two parame- 
ters ujq and U = A 2 /Mcuq to characterize basic properties 
of the system. The first parameter loq describes the adia- 
baticity. In the adiabatic limit of luq — > 0, the boson fields 
do not change in the imaginary time, that is, they behave 
as classical fields. In the opposite limit of ujq — » oo, the 
bosons react instantaneously to fermion motions. Be- 
tween these two limits, bosons with a finite u)q mediate 
a retarded effective interaction which is given by J7 G ff in 
Eq. (|l5|). The second parameter U describes the magni- 
tude of the effective interaction between fermions. Note 
that U = \U s(u> n = 0)| in this dispersionless case, since 
the bare impurity Green's function is given by the non- 
interacting one (|lg|). 

For a fixed value of luq, the system behaves quite dif- 
ferently in the regions with U <C 1 and U 3> 1 . For small 
values of U, fermions are nearly free and each lattice site 
is in an empty, a singly-occupied, or a doubly-occupied 
state with almost equal probability at half-filling (// = 0) . 
If we define the probability P(x) that the boson field x 
lies in the interval between x and x + Ax, P(x) shows a 
single broad peak centered at x — 0. Compared to this, 
if U becomes large, fermions strongly interact with each 
other, and a combined state between fermion and boson 
may be formed, which is called a small polaron. The 
polaron consists of double occupancy of fermions for the 
model with the interaction ( p"7| ) (bipolaron). Then, the 
probability P(x) displays a double peak at x = ±A/Mojq 
which corresponds to the doubly-occupied and empty 
states. Figure [l] shows this behavior by changing the 
value of U for the case of ujq = 0.5. The single peak of the 
probability P(x) appears for small U, while the double 
peaks are developed for U ^ 1 as shown in Fig. 0(a). At 
the same time, in Fig. |l|(b) , the probability of the double 
occupancy P-q increases from 1/4 for the noninteracting 
case to 1/2 for the situation in which the system consists 
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of only empty and doubly-occupied sites. The self-energy 
for fermion X is also enhanced by the effective interaction 
between fermions U. Figure ||(c) shows that the absolute 
value of the imaginary part of the self-energy as a func- 
tion of Matsubara frequency is strongly enhanced by U . 
Note that the data for uj n > 1/Ar contain no unbiased 
information. These clearly indicate the crossover from 
the weakly-correlated fermions in the sraaJl U region to 
the small polarons in the large U regionEO 

The similar crossover is found for other values of uiq. 
Figures || and || show the results for loq = 2 and 8, respec- 
tively. The value of U for the crossover, which we call U* 
hereafter, depends on the value of u)q. For example, for 
the case of u>o = 0.5 in Fig. |l|, the double-peak structure 
of P{x) appears at U ~ 1, on the other hand, it does 
not appear up to U ^3 for luq — 8. This can be under- 
stood as follows: In the limit of ujq — > oo, since the effec- 
tive interaction becomes spontaneous, U c s(t) = —U8(t^. 
the model maps onto an attractive Hubbard modeled 
in which the boson field coxEfisponds to the continuous 
Hubbard-Stratonovich field.Era In the Hubbard model, it 
is known that the continuous field develops a double-peak 
distribution at U ~ 3, which corresponds to the opening 
of the Hubbard gap in the case of a repulsive interaction.0 
On the other hand, in the opposite limit of loq — > 0, the 
effective interaction becomes constant in the imaginary 
time, U c s{t) = — U. This case is identical to an attrac- 
tive Falicov-Kimball model in the limit of a cpatinuous 
number of configurations for the static fields lirl In the 
adiabatic limit, fermions are localized at a smaller value 
of U since fluctuations of the boson field is smaller in this 
case than in the anti-adiabatic limit. Then, the splitting 
of the distribution of x should appear at a lower value of 
U. In the Falicov-Kimball model with a discrete-static 
field, the critical value of U is estimated to be l.Efl'Eil The 
finite value of luq can interpolate these two limits. Thus, 
the value of U* may change smoothly from U* ~ 3 in the 
limit of loq — » oo to U* ~ I in the limit of ujq — > 0. 



C. Dispersive Boson: Weak Coupling Limit 

Now we discuss the cases with a finite bosonic disper- 
sion; lui ^ 0. First, we study the weak coupling limit of 
W 3> We and U which has been studied by a perturbation 
theory.Q 

In this region, the finite width of the boson dispersion 
ui enhances the effective interaction between fermions. 
Figure ^(a) shows the bare impurity Green's function for 
boson T>q as a function of Matsubara frequency for var- 
ious values of lo\ for the case of uj = 0.5 and U = 0.16 
(A = 0.2). X>o is enhanced by the width of the dispersion 
which indicates that through the relation (15), the ef- 
fective interaction between fermions U e s is enhanced by 
This enhancement is also observed in the imaginary 
part of the fermion self-energy as shown in Fig. f|(b). At 
the same time, the probability of the double occupancy 



becomes large as shown in Fig. ||(c). These features 
are similar to those in Figs. |l[|3| when the parameter U 
increases in the small U region. These results can be 
understood using a perturbative argument in Sec. 



[V 



D. Dispersive Boson: Atomic Limit 

Next, we consider the limit of W <C ojq and U which 
has beep, .studied based on the so-called small-polaron 
theoryBo'El In this limit, the coherent band motion of 
fermions in Eq. (Q) is a perturbation on other terms of 
(||) and (^). The small-polaron theory is a perturbative 
approach from the atomic limit. The strong interaction 
between fermions and bosons leads to the formation of 
the small-po laron state as mentioned in the dispersionlcss 
case in Sec. Ill B . 

In this region, as the weak coupling case in Sec. [HI C , 



the effective interaction between fermions is enhanced by 
the finite width of the boson dispersion. Figure ||(a) plots 
the bare impurity Green's function for boson at zero Mat- 
subara frequency for loq — 8 and U = 9 (A = 24). A finite 
width of the boson dispersion W\ enhances V (uj n = 0). 
T>o shows the largest change at zero frequency as in Fig. 
||(a). At the same time, the absolute value of the imagi- 
nary part of the fermion self-energy increases as shown in 
Fig. ||(b). We plot here the data at the smallest Matsub- 
ara frequency to show the behavior clearly. The double- 
peak structure of the probability function P(x) shown in 
Fig. ||(a) at u>i = 0, which indicates the formation of 
the small-polaron state, does not change for u)\ within 
statistical errorbars. This suggests that the finite width 
of the boson dispersion enhances the effective interaction 
while the small-polaron state remains stable. These fea- 
tures will be discussed based on the small-polaron theory 
in Sec. Ifv. 



E. Dispersive Boson: Strong Fluctuation Regime 

Here we go beyond the perturbative regimes studied 
in Sec. IIIC and [II D. We consider the strong coupling 



case away from the anti-adiabatic limit, that is, U > W 
and loq ~ W. It is difficult to study this regime by any 
perturbative and analytical approach because of strong 
fluctuations. Our DMF method including the fluctuation 
effects is applied to this regime without any difficulty. 

Figure [} shows the results for luq — 0.5 and U = 2.56 
(A = 0.8). As shown in Fig. ^|(a), the absolute value 
of the bare impurity Green's function for boson T>q de- 
creases as u>i increases. The imaginary part of the self- 
energy for fermion also decreases its absolute value as 
shown in Fig. ^(b). At the same time, the probability of 
the double occupancy Pd decreases from 1/2 as shown in 
Fig. ||(c) . Figure ^|(d) shows that the double-peak struc- 
ture of the probability P(x) becomes unclear to merge 
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into a single peak. All these features exhibit that the ef- 
fective interaction between fermions U e s is weakened and 
the small-polaron state becomes unstable for u>\. This is 



a st riking contrast to the previous results in Sec. Ill C 



hid 



and 

havior in Sec. |TV 



We will discuss a physical picture for this be- 



In the intermediate region, we find a crossover as the 
value of uit increases. Figure ^ shows this crossover for 
u) = 0.5 and U = 0.64 (A = 0.4). For small values 
of u>i, we find a similar behavior as seen in Fig. ^; the 
bare impurity Green's function for boson is enhanced and 
both the absolute value of the self-energy and the dou- 
ble occupancy increase as uj\. However, for uj\ ^ 0.2, 
the behavior is reversed; all the three quantities turn to 
decrease as in Fig. ^. Therefore in this intermediate 
region, as the value of uj\ increases, the effective interac- 
tion between fermions is enhanced for small values of lo\ , 
however, turns to be weakened for large values of uj\. 

This crossover is closely related with complete soften- 
ing of the boson field. Figure || shows the effective fre- 
quency of the boson field uj* which is given by a pole of 
the Green's function for boson as 



uj* = a/ (uj - UJl) 2 + U(iuj n = 0), 



(24) 



where II is the self-energy for boson. The frequency uj* 
goes to zero at the value of ui\ where the crossover from 
the enhancement to the weakening of the effective inter- 
action exhibited in Fig. [?|. 



F. Phase Diagram 

We systematically investigate the crossover found in 
the previous section changing the parameter ujq and U. 
Figure || shows the values of uj* as a function of U for 
the cases of (a) wo = 0.5 and (b) loq = 2.0. For finite 
values of the width u>i, the frequency uj* goes to zero in 
both cases. We determine the crossover values of U by 
this complete softening of uj* for ujq and uj\. 

Figure [l0| summarizes the phase diagram for the 
crossovers determined by the above criterion. This indi- 
cates the boundary between the weak-fluctuation and the 
strong-fluctuation regimes as discussed in Sec. [V. The 



most important point in this phase diagram is that espe- 
cially for large uj^, the energy scale of U for th is cro ssover 
is quite different from U* determined in Sec. IIIB. This 
suggests that there is another parameter which controls 
the onset of strong fluctuations as discussed in the next 
section IV. 



IV. DISCUSSION 

In this section, we discuss the results obtained in Sec. 
[II. Pcrturbative arguments are applied to discuss the en- 
hancement of the effcfitivc interaction between fermions 
in the weak-couplinga and the atomic regions Ja^El In 



the strong coupling regime away from the anti-adiabatic 
limit, the weakening of the effective interaction and the 
instability of the small-polaron state are discussed as a 
consequence of the strong fluctuations of the boson fields 
accompanied by complete softening. The phase diagram 
is examined to clarify the parameters which control the 
onset of the strong fluctuations. 

In the weak coupling region, the dispersion width uj\ 
enhances the effective interaction between fermions in our 
DMF solutions. The absolute value of the self-energy for 
fermion is enhanced. A first-order perturbation in the 
coupling parameter A concludes that the self-energy £ be- 
comesJarger as the width uj\ increases since 2?o increases 
as uj\a Thus, perturbation theory suggests that the bo- 
son dispersion increases the effective interaction between 
fermions in the weak coupling limit. This enhancement 
can be understood intuitively as follows: In the weak cou- 
pling region, the density of states for both fermion and 
boson are not altered drastically by the fermion-boson 
interaction; a rigid-band picture should be justified. For 
a finite uji , the band edge of the boson density of states 
is lowered linearly. Then, the effective interaction (15) is 



mainly mediated by the bosons near the band edge as 



U eS {iUJ n ) 



A 2 



(iuJn) 2 - (LOO - ^l) 2 



(25) 



Thus the absolute value \U e g\ becomes large as the value 
of uji increases. Therefore, the enhancement of the effec- 
tive interaction in the weak coupling region can be un- 
derstood as the dec rease of the effective boson frequency. 
Our results in Sec. [II C are consistent with this pertur- 
bative argument. 

We now turn to the atomic limit, W <C ujq and U . Now 
the fermion hopping term in Eq. (|7|) is a perturbation to 
the terms (||) and (g) . If we apply the canonical transfor- 
mation to diagonalize the imnerturbed terms according 
to the small-polaron theoryBQu we obtain the expression 
of the Hamiltonian as 

q ^ ' ia 

+ T,^ c lc 3a X}X v (26) 

ij,a 

where A is the stabilization energy of polarons given by 

(27) 

and the operator Xi takes the form 



A = E-> 

q ^ 



Xi = exp 



A 



(28) 



The third term of the Hamiltonian ( p6[ ) indicates that the 
hopping occurs not as a bare fermion but as a combined 



6 



object between fermion and boson. Each fermion is asso- 
ciated by a local boson. This is called the small-polaron 
state. 

In the unperturbed state (ty = 0), the polarons are 
almost localized in real space with the stabilization en- 
ergy A given by Eq. (^). When one increases the width 
of the boson dispersion u>i , the stabilization energy A in- 
creases. This makes the polarons more strongly localized. 
The derealization of the polarons is a second-order per- 
turbation in i-term in Eq. ( p6| ) since the polarons consist 
of the double occupancy of fermions in this model. The 
hopping matrix by the second-order process is suppressed 
by u>i because the intermediate state in the perturbation 
costs the energy A. The operator X t does not change the 
result in this limit of ujq 3> W. Therefore a finite width of 
the boson dispersion suppresses the band motion of the 
polarons and enhances the effec tive interaction in this 
limit. Our results in Sec. [II D| are in good agreement 
with this argument based on the small-polaron theory. 

We now discuss the strong coupling reg ime a way from 
the anti-adiabatic limit studied in Sec. HIE. In this 



regime, the polaron state is formed by the strong cou- 
pling, however the boson fields are loosely bound to the 
fermions due to a finite u>o, compared to the previous 
atomic limit with ujq 3> W where the boson fields re- 
act instantaneously to fermion motions. This leads to 
large fluctuations in the boson fields by the hopping of 
fermions. Thus this regime is characterized by these 
strong fluctuations, which is the reason why any per- 
turbation cannot be applied. These fluctuations may in 
turn accelerate the derealization of fermions through the 
mutual feedback effects of the fermion-boson coupling. 

A finite width of the boson dispersion uj\ introduced 
in our calculations increases the fluctuations of the boson 
fields. The bosons are not localized and gain their kinetic 
energy through the dispersion. By tuning the width wi, 
we can control the fluctuat ions of the boson fields by 
hand. Our results in Sec. Ill E clearly exhibited that 



the effective interaction between fermions is weakened by 
u>i . This is considered to be a consequence of the strong 
fluctuations of the boson fields enhanced by u>i which 
tend to make fermions more delocalized. This behavior 
is elucidated for the first time by our method which fully 
includes the mutual feedback in many-body systems. 

In the intermediate coupling region, a sharp crossover 
was found by changing the value of uj\ in Sec 



HIE 



For small ui, the effective interaction is enhanced by 
uo\. Since the fluctuations are small there, this may 
be smoothly connected to the behavior discussed in the 
weak-coupling or the atomic regions. The effective fre- 
quency ui* defined by Eq. (^J) becomes small but re- 
mains finite as in the perturbative regime although the 
reduction of u* is large and nonlinear in this nonper- 
turbative regime. When the value of lo± becomes large 
enough to soften the boson field completely (u>* — > 0), 
fluctuations play a crucial role to enhance the derealiza- 
tion of fermions. The boundaries in the phase diagram in 
Fig. nn are the crossovers between the weak-fluctuation 



and the strong-fluctuation regimes. 

In the dispersionless case in Sec. Ill B , we have found 
another crossover by the formation of the small polaron 
which is, for instance, characterized by the development 
of the double-peak structure in the probability P{x). The 
critical value of U for this crossover, U* , changes from 
U* ~ 3 in the limit of ujq ^> W (anti-adiabatic limit) to 
U* ~ I in the limit of ojq <C W (adiabatic limit) as shown 
in the schematic phase diagram in the plane (1/ujq,U) 
in Fig. |ll] (dotted gray line). On the other hand, the 
crossover to the strong fluctuation regime in Fig. |l^ ap- 
pears at much larger values of U than U* especially in 
the anti-adiabatic regime with large but finite ojq. This 
strongly suggests the importance of another energy scale 
to characterize the strong fluctuation regime which is not 
clearly found in the dispersionless case. 

The importance of such characteristic energy scale has. 
been pointed out also in a previous mean-field study.E£l 
The parameter is defined by the ratio of the fermion- 
boson interaction to the spring constant of boson fields, 
r/ = A/ MuJq . In the case of rj < 1 , since the fermion-boson 
interaction is weak compared to the stored energy in the 
boson field, the single-boson process should be impor- 
tant. In the case of r\ > 1, the fermion-boson interaction 
is strong enough to excite a large numbers of bosons, i.e., 
multibeson processes become important. In the previous 
study,c3 the importance of fluctuations of the boson fields 
has been suggested in the latter multiboson regime. 

In the plane (l/t<;o, U), the crossover between the 
single-boson and the multiboson regimes occurs at rj = 1 , 
i.e., for U = luq and is shown in Fig. [ll] (solid line). The 
line of 1] — 1 becomes much larger than U* in the anti- 
adiabatic regime. If we plot these values of U(j] = 1) 
on the axis of u>i = in Fig. 10, the crossover bound- 
aries seem to be smoothly connected to these values in 
the anti-adiabatic region. We demonstrate this behavior 
in Fig. 12 for luq — 4 and 2 (gray lines). This indicates 
that in the anti-adiabatic regime, the line of rj = 1 cor- 
responds to the crossover between the weak- fluctuation 
and the strong fluctuation regimes in the dispersionless 
case. 

On the other hand, in the adiabatic regime with small 
but finite u>q, the value of U for the boundary in Fig. [l(] 
is not smoothly connected to the value of U for 77 = 1. 
For instance, in the case of uj — 0.5, 77 = 1 gives U =1/4 
which is much smaller than the boundary. This suggests 
that the condition of 77 > 1 does not characterize the 
strong fluctuation regime in this adiabatic region. 

To understand this behavior in the adiabatic regime, 
let us discuss the adiabatic limit of ujq — > 0. In this limit, 
the boson fields behave as classical fields which do not 
fluctuate in the imaginary time direction. Boson fluctu- 
ations only come from the fluctuation of the value of the 
field x which is constant in time. Then even if the system 
is in the region of 77 > 1, the fluctuations of the boson 
fields are small when U is smaller than U* since the fields 
feel a deep single- well potential as indicated in the prob- 
ability P(x). The boson fields begin to fluctuate when 
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U becomes comparable to U* where the potential for x 
softens around x = 0. Therefore in the adiabatic limit, 
the strong fluctuation regime should be characterized not 
by 77 > 1 but by U > U* . 

Based on this argument, if we plot the value of U* 
on the axis of u)\ — 0, the crossover boundaries in the 
adiabatic regime seem to be smoothly connected to these 
values for the cases of loq — 0.5 and 0.25 in Fig. [l2| 
(dotted lines). 

These results reveal that the boson fluctuations which 
are strong enough to delocalize fcrmions appear in a 
different way in the anti-adiabatic and the adiabatic 
regimes. In the anti-adiabatic regime (coq > W), the 
small-polaron state is formed at U is U* . In the region 
with U > U* and rj < 1, however, the fluctuations of the 
boson fields are small in the sense that the single-boson 
process contributes mainly and that the effective boson 
frequency is finite. If rj becomes larger than 1, the boson 
field is softened and the boson fluctuations play a crucial 
role through the multiboson process. On the other hand, 
in the adiabatic regime (lo < W), the fluctuations do 
not become large until U ~ U* even if rj is larger than 1. 
The fluctuations are mainly the classical origin there. 

To summarize, the strong fluctuations of the bo- 
son fields become important only when the conditions 
U > U* and rj > 1 are both satisfied. These condi- 
tions are shown as the hatched area in Fig. [n]. This 
area is strongly-correlated region for both fermions and 
bosons. The criterion for the formation of small polarons, 
U ~ U* , corresponds to a competition between the ki- 
netic energy of fermions W and the effective interaction 
U. The line of U* may be modified according to a spe- 
cific form of the fermion-boson coupling. On the other 
hand, the criterion rj ~ 1 corresponds to a competition 
between the stored energy of the boson field and the cou- 
pling energy to fermions. Thus the hatched area in Fig. 
|l l| is the region where correlations become strong in both 
standpoints of fermions and bosons. 

A subtle problem remains open about the boson soft- 
ening in the dispersionless case. As shown in Fig. ||, 
when ui is zero, the effective frequency to* becomes very 
small but remains finite for large values of U . We note 
that there are finite-temperature effects; lu* is suppressed 
more strongly for lower temperatures. Unfortunately we 
cannot conclude in this study whether lu* goes to zero 
even when lu± = 0. In this dispersionless case, the boson 
density of states is a delta function at e = u)q, which is 
special since the shape of the boson density o f stat es at 
the bottom is important as mentioned in Sec. HI A. For 
instance, the step-like singularity in the two-dimensional 
density of states might prevent the boson field from com- 
plete softening at finite temperatures. Though further 
studies are necessary for the property of DMF equations 
for various types of density of states, we believe from the 
results in Figs. |ll| and [l^ that the nonlinear suppression 
of lu* is relevant to strong fluctuations of bosons and that 
there are two important energy scales even when uj± = 0. 



V. SUMMARY AND CONCLUDING REMARKS 

We have investigated the effects of the boson disper- 
sion in a system of dynamical mean-field equations de- 
scribing coupled fermion-boson systems. The analysis of 
the equations revealed that the boson dispersion plays 
a crucial role in a wide region of parameters. By intro- 
ducing a parameter for the width of the dispersion in the 
model, we can control the fluctuations of the boson fields. 
To handle the boson fluctuations and the feedback ef- 
fects, we have extended the dynamical mean-field theory 
to determine the Green's functions for both fermion and 
boson in the self-consistent way. In the ordinary frame- 
work for the dispersionless case, the channel for the boson 
Green's function is frozen in the sense that the bare im- 
purity Green's function is fixed and unrenormalized from 
the noninteracting one. The renormalization of the bare 
impurity Green's function for boson is very important 
since the bare impurity Green's function is directly re- 
lated to the effective interaction between fermions. The 
equations in the extended dynamical mean-field theory 
are solved by using quantum Monte Carlo technique. 

The main result in the models with dispersive bosons 
is that in the strong coupling regime away from the anti- 
adiabatic limit, the fluctuations of the boson fields be- 
come relevant to accelerate the derealization of fermions. 
The effective interaction between fermions is weakened as 
the width of the boson dispersion increases in this regime. 
This behavior is explicitly shown for the first time by our 
method which fully includes the mutual feedback effects. 
The crossover to this nonperturbative regime is closely 
correlated with softening of the boson field. We have ex- 
amined the phase diagram where this strong fluctuation 
occurs by tuning the coupling parameter and the width 
of the dispersion. The strong fluctuations to delocalize 
fermions become relevant when the small-polaron state is 
formed and the multiboson processes become important. 
The small polarons become stable when the effective in- 
teraction between fermions overcomes the fermion band 
energy. The multiboson regime is characterized by a cou- 
pling parameter larger than the boson energy. Thus the 
strong fluctuation regime is the strong correlated region 
for both fermions and bosons. As the coupling parameter 
increases, the boson fluctuations appear in a different way 
between in the adiabatic and the anti-adiabatic regimes. 
In the adiabatic regime, the fluctuations are mainly clas- 
sical which are enhanced by the softening of the potential 
for the boson fields in the formation of the small-polaron 
state. On the other hand, in the anti-adiabatic regime, 
the small polarons are formed in the single-boson regime, 
where the dynamical fluctuations are small and the ef- 
fective boson frequency is finite. The boson fluctuations 
do not play a crucial role until the system enters in the 
multiboson regime by complete softening of bosons. 

The onset of the strong fluctuations occurs near the re- 
gion where the boson degrees of freedom soften. In this 
paper we have studied the DMF equations in the absence 



8 



of any freezing of the boson degrees of freedom. These 
effects together with generalizations to states with differ- 
ent symmetries and other generalizations are currently 
under investigation. 

Our results imply that the behavior of the boson fluc- 
tuations may depend on the specific form of the bo- 
son density of states. Different forms of the density of 
states should be tested in the present DMF framework 
in a future study. Especially we are interested in the 
two-dimensional case with a step-like singularity at the 
edge which might be free from complete softening at fi- 
nite temperatures. Boson fluctuations in this case may 
lead to light-mass bipolaronic states, which would give 
some insights into the high-temperature superconductiv- 
ity in Cu-oxide materials where fermions strongly couple 
with spin fluctuations. We plan to understand this two- 
dimensional case in a later publication. 

The dynamical mean-field equations allow us to vary 
the width of the boson dispersion in the calculations. 
This reveals the interesting properties in the strong fluc- 
tuation regime. The tuning of the electronic bandwidth 
has been the subject of a great deal of theoretical and ex- 
perimental workB3. Our work suggests the possible inter- 
est of varying the boson dispersion experimentally even 
though this may be easier in systems where the bosons 
are spin fluctuations whose dispersion determined by ex- 
change interactions can be controlled more easily than 
optical phonon dispersions. Another possibility may be 
the realization of the dynamical mean-field theory in a 
random model. 

There are many materials which satisfy the above con- 
ditions for the strong fluctuation regime. In many physi- 
cal situations, the fermion bandwidth W is large or com- 
parable to ujq, which makes possible to access to the 
strong fluctuation regime by a relatively weak coupling. 
Our method provides a powerful theoretical tool to exam- 
ine the physical properties in this regime. We can apply 
it to more realistic models including orbital degrees of 
freedom of electrons, different normal modes of phonons, 
or interactions between fermions. Such extensions are 
now under investigation. 
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FIG. 1. Results for the dispersionless model with loo = 0.5 
at (3 = 8; (a) the probability function of the boson fields 
x, (b) the probability of the double occupancy, and (c) the 
imaginary part of the self-energy for fermion as a function of 
Matsubara frequency. In (a), the typical errorbars are shown 
at the peaks of the distributions. The lines in (b) and (c) are 
guides to eye. 
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FIG. 2. Results for the dispersionless model with u>o = 2. 
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FIG. 10. The phase diagram of the crossover for various 
values of ljq. The boundaries are determined by the complete 
softening of the effective boson frequency. The lines are guides 
to eye. 
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FIG. 11. Schematic phase diagram for the dispersionless 
case. The critical value for the formation of the small-polaron 
state U* is shown as the dotted gray line. The solid line 
indicates the boundary r) = 1. The hatched area is the strong 
fluctuation region. See the text for details. 
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FIG. 12. The values of U and uj\ at which the crossover 
from the weak-fluctuation to the strong-fluctuation regime 
takes place. Notice that these values extrapolate smoothly 
to the values of U at r\ — 1 for ujq = 4 and 2 (anti-adiabatic 
regime) and to the U* for luq — 0.5 and 0.25 (adiabatic regime) 
in the limit of uj\ = 0. See the text for details. 
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